Spin-Orbit mediated spin relaxation in monolayer M0S2 
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We study the intra-valley spin-orbit mediated spin relaxation in monolayers of M0S2 witliin a two 
bands effective Hamiltonian. The intrinsic spin splitting of the valence band as well as a Rashba-like 
coupling due to the breaking of the out-of-plane inversion symmetry are considered. We show that, 
in the hole doped regime, the out-of-plane spin relaxation is not very efficient since the spin splitting 
of the valence band tends to stabilize the spin polarization in this direction. We obtain spin lifetimes 
larger than nanoseconds, in agreement with recent valley polarization experiments. 



I. INTRODUCTION 

Among all the quasi two dimensional crystals that 
have become popular since the appearance of graphenCj^^ 
monolayer unit cells of molybdenum disulfide^- (M0S2) 
and other dichalcogenides are particularly attractive be- 
cause the existence of an electronic gap which makes 
those systems an excellent candidate for nanoelectronics 
devices.'^''^ Whereas multilayer samples presents an indi- 
rect band gap, monolayer M0S2 has a direct gap across 
the inequivalent K and K' points of the hexagonal Bril- 
loin zone (BZ), which makes this material of particu- 
lar interest for optoelectronic applications. Of special 
interest is the promising applications of M0S2 for spin- 
tronics and valleytronics devices .l^ltiS! Jn fact, the strong 
spin-orbit (SO) interaction, together with the absence of 
inversion symmetry in monolayer samples, splits the va- 
lence bands by ~ 150 meV in two spin flavors. This 
splitting is essential for spintronics and optoelectronics 
applications and, importantly, it has different signs at 
each valley due to time reversal symmetry. This fact 
allows to control the valley population by optically excit- 
ing the monolayer samples with circularly polarized light, 
as it has been demonstrated experimentally,'^ what is 
generically called spin-valley couplingii - 

A crucial role for the efficiency of the valley polariza- 
tion is played by the spin lifetime Tg of the system, which 
must be longer than ^ 10 ns for realistic applications. 
Interestingly, coherence times of this order have been ex- 
perimentally measured for single layers of MoS2.^ There- 
fore, understanding the spin relaxation mechanisms of 
this material is essential. Due to the strong SO coup ling 
and the presence of disorder, both Elliot- YafePE^l and 
D'yakonov-PerelfSEH mechanisms are expected to play 
a role. Here we perform a systematic calculation of the 
intra-valley SO mediated spin relaxa tion rates by means 
of the Mori-Kawasaki formula which is the appro- 
priate framework to treat both mechanisms on the same 
footing. Our results show that the intrinsic reflection 
symmetry of the system with respect to the out-of-plane 
direction in combination with the large spin splitting of 
the valence band allows spin lifetimes for the out-of-plane 
polarization larger than nanoseconds, in agreement with 
the experiments.^ 

The manuscript is organized as follows. First, we 
present the two bands effective model that we employ 
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FIG. 1: a) Top view of the lattice in real space. Red and 
blue dots represent the Mo and S atoms, respectively, b) 
Brillouin zone, with the corresponding reciprocal lattice wave- 
vectors. The inequivalent -|-K and — K points are shown, 
c) Conduction and valence bands at ±K points within the 
effective model discussed in the text. The SO splitting is not 
considered in this figure. The numbers inside the kets which 
label the bands express the phases picked up the Bloch wave 
function under symmetry operations of the crystal structure, 
see Tabs. H] and ini 



to perform the calculation, including an exhaustive dis- 
cussion on the microscopic origin of the intrinsic and 



Rashba-like SO couplings. In Sec. Ill we compute both 
the in-plane and out-of-plane spin relaxation rates associ- 
ated to these couplings. Finally, in Sec |IV| we summarize 
our results. 



II. THE MODEL 

A. Two bands effective model 

The transition metal dichalcogenide M0S2 is com- 
posed, in its bulk configuration, of two-dimensional S- 
Mo-S layers stacked on top of eachother, coupled by week 
van der Waals forces. The Mo atoms are ordered in a 
triangular lattice, each of them bonded to six S atoms 
located in the top and bottom layers, forming a sand- 
wiched material. A top view of the lattice is shown in 
Fig. [l]a). Like in graphene, the weak interlayer coupling 
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Csh = C3 X ah 
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TABLE I: Character table of Czh- We denote = e 3 . 
Notation from Ref. [20l 

makes possible to exfoliate this material down to a single- 
layer!^ The electronic band structure of M0S2 changes 
from an indirect band gap for multilayer samples, to a 
direct gap semiconductor for single-layers, located at the 
two inequivalent ±K points of the BZ [see Fig. [l]b)]P 



Irreps 


C3 


Oh 


Mo 


s 


A! 


1 


1 


75 (d^2_y2 ± id^y) , 
75 [Px T iPy) 


^ {Px ± iPy) (S) 


A!' 


1 


-1 


1"^ 


75 (Px ± iPy) (as) 


E' 




1 




75 (P:c =F iPj,) (s) 


E' 




1 




Pz (as), s (as) 


E" 




-1 


Pz 


75 {Px =F iPa) (as) 


E" 




-1 




Pz (s), s (s) 



TABLE II: Classification of the Bloch wave functions at BZ 
corners according to the 6 irreducible representations of Csh- 
The sign ± corresponds to ±K points. In the case of S 
atoms, both symmetric (s) and anti-symmetric (as) combi- 
nations with respect to Oh of the orbitals of the top and the 
bottom atoms are considered. The second and third column 
contain the phases picked up by the wave function when a 
2tt/3 rotation or a mirror reflection is performed. 



Being interested on single-layer samples, we compute 
here the spin relaxation rates within the 2 bands effec- 
tive k • p Hamiltonian at the two inequivalent corners of 
the BZ. The hexagonal D^ii symmetry of the monolayer 
crystal determines the explicit form of the Hamiltonian. 
Note that the D^h point group can be expressed as the di- 
rect product of D3, which contains the identity, the two 
counterclockwise rotations by 27r/3 and 47r/3, and the 
reflections across the three in-plane axis which connect 
Mo and S atoms in a top view of the lattice [see Fig. [l] 
a)], and cr/j, which contains the identity and the inversion 
across the out-of-plane axis. The group of the wavevec- 
tor at ±K is C3/1 = C3 x cr/j (the reflections across the 
in-plane axis swap the two inequivalent K points). The 
character table of this group can be found in Tab.|l) Note 
that there is no two-dimensional irreducible representa- 
tions of Csh, so band touching at Dirac points are not 
protected by symmetry. Each band at ±K can be la- 
beled by two quantum numbers which express the phase 
picked up by the corresponding Bloch wave function at 
these points under a rotation by 27r/3 and a reflection 
across the out-of-plane axis respectively, or equivalently, 
the irreducible representation of C3/1 associated to the 
Bloch wave function. These phases depend, of course, on 
the orbital character of these states. In Tab. Illla classifi- 
action of the Bloch wave functions at the BZ corners for 
Mo orbitals (s, p, d) and S orbitals (s, p) is given. 

Around ±K points the conduction and valence bands 
are mostly made of d orbitals coming from Mo. In partic- 
ular, the orbital weight of the conduction band is essen- 
tially d^z^-r^, which belongs to the E' irreducible repre- 
sentation of C3/1, whereas in the case of the valence band 
is mostly the real combination of d^2_y2 and dxy which 
belongs to A' . We restrict our analysis to a minimal low 
energy model that accounts for the above configuration. 
Up to first order in k, the effective Hamiltonian readsf^ 

Ho = atiraxka: + ayky) + (1) 



where the Pauli matrices at operate in a space of 2- 
component Bloch functions — {ip\w^ ■ 
Here |w'^(l),-fl(-|-l)) labels the symmetry properties of 
the conduction (valence) band state, and t = ±1 corre- 
sponds to valley ztK.l^Hl 



In order to discuss the effective SO coupling within 
this model we must introduce Pauli matrices associated 
to the spin degree of freedom. Importantly, we are now 
introducing a pseudovector in the 3-dimensional space, 
meaning that the operators which contain Sz are even un- 
der ah whereas the operators which contain the in-plane 
component of the spin are odd. Unless the ah symmetry 
is expressly broken, our effective model may only contain 
terms with Sz- Therefore, the ah symmetry of the crys- 
tal structure protects the out-of-plane spin component. 
Then, the intrinsic SO coupling terms read in general: 

nL'i^XcT^Sz + X.r^^^Sz (2) 

where I is the identity matrix acting in the space of 2- 
component Bloch functions. The absence of a center of 
inversion in the crystal implies the spin splitting of the 
bands. Here Ac and A^, are the splittings of the con- 
duction and valence bands respectively. Although both 
splittings are allowed by symmetry, it is important to 
notice their different microscopic origins due to the dif- 
ferent orbital character of the bands. We consider an 
intra-atomic SO Hamiltonian of the form A50L • s for 
the d orbitals of Mo, where L is the angular momentum 
operator. As we show schematically in Fig. [2j the spin 
splitting of the valence band is the result of a first order 
process. However, the splitting of the conduction band 
is associated to second order processes which involve vir- 
tual transitions into states wich belong to the A" and E" 
irreducible representations. Therefore, 

Aint = At, 3> Ac, (3) 
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FIG. 2: Sketch of the microscopic processes which lead to 
the effective SO coupling terms discussed in the text, a) First 
order processes which lead to the splitting of the valence band, 
b) Second order processes associated to the splitting of the 
conduction band, c) Second order processes which lead to a 
Bychkov-Rashba coupling when an symmetry is broken. 



and we neglect from here on the SO splitting of the con- 
duction band. 

The values of the model parameters can be ex- 
tracted from experiments as well as from first principle 
calculations;^^'^^ Here we take at — 3.51 eVA, A = 1.66 
eV and 2Amt — 0.15 eV, with the lattice constant a = 
3.193 A. 



B. Rashba effect 

If the CT/i symmetry is broken then is not longer a 
good quantum number. In that situation a coupling with 
the in-plane components is possible. A Bychkov-Rashba 
coupling^^ appears, which in the context of the two bands 
effective model reads: 



(4) 



This coupling is the result of second order processes as 
in the case of the spin splitting of the conduction band, 
although first order in the SO interaction as it is shown 
in Fig.|c). 

The Rashba effect in 2D systems is usually attributed 
to dipolar transitions induced by the aplication of an elec- 
tric field in the out-of-plane direction. In our case, such 
field would induce dipolar transitions between the Mo 
d3z2_^2 orbitals of the conduction band and Mo or- 
bitals of bands at much higher energies. Then, the SO 
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FIG. 3: Spin susceptibility x as a function of the chemical 
potential /i, for different values of F, calculated from Eq. ( 15 1 . 



Solid black line corresponds to F = 0.001 eV, blue dashed line 
to F = 0.01 eV, and red dotted line to F = 0.1 eV. 



interaction would induce transitions between these states 
and the valence band flipping the spin. However, the or- 
bital weight of the valence band in Mo p orbitals is very 
small, so this coupling is expected to be very weak. Nev- 
ertheless, this picture changes if orbitals from S atoms 
are also taken into account. For instance, if we consider 
the application of a gate voltage Vgate, which is necessary 
in order to induce charge carriers in this system, then we 
would have different on-site energies for the p orbitals 
of the top and bottom S atoms. This turns into a non 
zero hybridization between valence band states 1 1,-1-1) 
and states |1, —1) of higher bands proportional to Vgate- 
Then, the SO interaction induce transitions between d^.^ , 
dy2 orbitals of these bands and d^^^^^^ of the conduction 
band, flipping the spin. This kind of processes is the 
one depicted in Fig. [2]c). Since the orbital weigth of S 
p orbitals in these bands is less than the 20 %, we can 
estimate an upper limit for this coupling of the form: 
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(5) 



where e|i,-i) represents the energy of the |1,— 1) {A") 
band involved in the calculation. 



C. Disorder and Mori-Kawasaki formula 

We compute the spin relaxation rates by using 
the Mori-Kawasaki formula. Originally, the Mori- 
Kawasaki formula was deduced in order to compute the 
broadening of the signal peak in an electron spin res- 
onance experiment due to the breaking of the SU{2) 
spin symmetry of the system. As it has been shown 
recently,'^ it can be related with the inverse of the spin 
lifetime — H/ts. The Mori-Kawasaki formula treats 
the SO coupling terms as perturbations to the electronic 
Hamiltonian, something that in principle is valid since 
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in our model Xint, Xext are at least one order of mag- 
nitude smaller than the other energy parameters of the 
model, namely the gap A and the bandwidth 2t. The 
spin lifetime is computed as Ts- 



-— lini Im "^"^"^^ ^ — ^ 



(6) 



where x is the spin susceptibility and Xaa^ i'-^) the 
Fourier transform of the response function: 



(7) 



with A = [n^^,s+] and s± = (s^ ± isy)/2. The spin 
susceptibility is defined as: 



dH 



(8) 



H=0 



where H is the field strength associated to a Zeeman term 
in the Hamiltonian "Hz = —g^sHsz.^'^ The expectation 
values in Eqs. ([t]) and ([s]) are referred to the Hamiltonian 
without the SO coupling terms. 

For SO mediated spi n rel axation mechanisms, par- 
ticularly th e Elli ot- Yafeij^^EH ^nd the D'yakonov-Perel' 
mechanisms,'i^'i^ a relevant parameter of the theory is 
the amount of disorder F = h/r-p, where Tp is the lifetime 
of the quasiparticles with momentum p. Disorder is in- 
troduced in our model in a phenomenological manner, by 
adding the imaginary self-energy iF/2 to the Matsubara 
Green's operator associated to the free Hamiltonian Q: 



G(k,zu;) = ^ ^ G„(k,*c.) 

a=±l 

where we define: 



I + a \ afk ■ 3 + i^'^z^ /^k 
(9) 



be calculated in terms of the spectral functions ^^(k, w), 
defined from the retarded version of the Green's functions 



of Eq. ( 10 1 as 



(k,w) = -2 ImGf(cj,k) = 
Then, we can write: 



(12) 



/oo 
— ni^(a;)Aa(k,a; + sg^siJ). (13) 
-oo ^''^ 



a=±l 



where np{uj) is the Fermi-Dirac distribution function and 
the factor 2 accounts for the valley degeneracy. From the 
definition of Eq. (Is]) we get: 



1 



^ N ^ 

Q=±l k ' 



^Aa(k,C.)(- 



dnp{ijj) 



(14) 



In the zero temperature limit (T <^ Tp, where Tp is 
the Fermi temperature) we can approximate — ^"gj"'* ~ 
S{lj). The sum in k can be written as an integral through 
the standard procedure ^ — )■ / d^k, where Ac is 
the area of the unit cell. The isotropy of the dispersion re- 
lation allows to integrate in angles straightforwardly and 
to write down the remaining integral in |k| as an integral 
in energies. At this point, it is necessary to introduce 
an energy cutoff D for the effective model, which can be 
related with the area of the unit cell as D — at^i: j A^,- 
After some algebra the spin susceptibility can be written 
as: 



7r2a2t2 2^ 

a=±l 



dx- 



A 

2D 



4Z52 



(15) 



Ga(k, iw) 



ack + /U + ir/2 



(10) 



Here ek = \/ a^t^jkp -|- A2/4 is the dispersion relation 
of conduction (a. = +1) and valence (a = —1) bands 
of our effective model and /z is the chemical potential. 
Note that the valley index r is omitted since we do not 
consider short-range scatterers which could connect both 
valleys. Therefore, inter-valley relaxation processes are 
beyond the scope of this work. We also neglect the mo- 
mentum dependence of F, so it enters just as a parame- 
ter which in principle can be determined from transport 
experiments.'^ 

The expectation value of the z-component of spin is: 



(11) 



where TV is the number of unit cells and is the occu- 
pation number of quasiparticles with momentum k and 
spin s in the presence of the Zeeman term Hz- This can 



The spin susceptibilty as a function of the chemical po- 
tential is shown, for different values of F, in Fig. [3j If 
we drop logarithmic terms that appear in Eq. (15]) after 
integration in x, we obtain a simple analytica 
for X which is valid for ^ ^ D: 



formula 



^ 2fiAc \ - 

a=±l 



arctan 



2fi ~ aA 



(16) 



We compute now the numerator of Eq. ([6|. The cal- 
culation is easily performed in the Matsubara frequency 
domain. As we said before, we neglect the momentum 
dependence of F. Thus, vertex corrections in the compu- 
tation of x_4_4i {itL)) can be ignored. Then, we can write: 



k iv a.. a' 



(17) 

where /3 is the usual thermal factor and faa' (k) is defined 
as: 
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faa' (k) = -Tr 



A- \ I + a- 



A , 



, atk. ■ a 



(18) 



The trace is performed in the space of 2-components 
Bloch functions, and the vahey degeneracy has already 
been taken into account in this definition. The sum in fre- 
quencies can be performed easly by using the Lehmann 
representation in terms of the spectral functions intro- 
duced before. After the sumation and the analytical con- 
tinuation we have for the imaginary part of Xaa'^ i^)- 



de 
47r 



Aa (k, e + uj) Aa' (k, e) [np (e) - np {e + w)] 



(19) 



Hence, in the a; — >■ limit we obtain: 



k a, a' 



f°° de 
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(20) 



After the same approximations as before we can write, 
in the zero temperature limit: 



T, = 



1 A 



Attx (27r) 



^ E / '^^^f^-' (^) (1^' 0) (k, 0) 



(21) 

The remaining part of the paper will be devoted to the 
estimation, using Eq. (21), of the spin relaxation in the 



different scenarios which are relevant for M0S2. 



III. SPIN RELAXATION 



A. In-plane spin relaxation 
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FIG. 4: In-plane spin lifetimes as a function of the carrier 
concentration. Top: Electron doping. Bottom: Hole doping. 
Dashed black line corresponds to F = 0.001 eV, dotted blue 
to F = 0.01 eV, and solid red F = 0.1 eV. Inset: In-plane 
spin lifetimes for electron concentrations in double logarith- 
mic scale. Notice the different time scale in the top and bot- 
tom panels. 



We start by computing the in-plane spin relaxation 
rate due to the intrinsic SO coupling. For this, we use 
Eq. (21 ) with A = Xmt [I — o'z) s+, which leads to: 



U^, (k) = 2AL !-(« + «') 



2ek 



aa'^l (22) 



Then, one sees that the in-plane relaxation rate can be 
written as the sum of two contributions, one coming 
from intra-band transitions and the other from inter- 
band transitions: 



(23) 



where: 

-^intra 



E _ 

Q=±l " 2D 



dx 



dx 



{x-y) +^ {x + ^l) + 



4 



{x - a^i) + i^J 1^(2; - an) + ^ 

(24) 



The intra-band transitions account the D'yakonov- 
Perel' processes, whereas the inter-band term leads to 
the Elliot- Yafet contribution. This is more clear in the 
doped regime |/.t| > A/2. If we drop logarithmic correc- 
tions in the above integrals, as we did in order to get 
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Eq. (16), we arrive at: 



\2 

-^intra ^ int 
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(25) 



The inter-band transitions lead to an Elliot- Yafet con- 
tribution characterized by the linear scaling between the 
spin lifetime and momentum scattering time T™*'^'' oc 
r. The intra-band transitions, however, leads to the 
D'yakonov-Perel' mechanism, characterized by F™*'"" c>c 
r~^. This mechanism is clearly the dominant one, as ex- 
pected from symmetry considerations, due to the absence 
of a center of inversion in the crystal structure. Assuming 
that r <C /i, A we have: 



^intra 

in 
-winter 
' in 



1 - 
1 + 



2fj, 



r/2 



(26) 



and therefore r™*''"/r™*'^'' ^ 1 unless the chemical po- 
tential lies at the bottom of the conduction band. It is 
important to note that the D'yakonov-Perel' mechanism 
is clearly electron-hole asymmetric due to the different 
spin splittings of the conduction and valence bands. 

These features are clearly shown in Fig. |4j where the 
in-plane spin lifetime is computed numerically. We see 
that the D'yakonov-Perel' mechanism is clearly domi- 
nant for hole dopings. From mobilities reported in trans- 
port experiment^- we deduce F « 0.02 eV, and therefore 
Tin ~ '^f'-^/^int ~ 5 fs. For electron concentrations, it is 
interesting to note the crossover from D'yakonov-Perel' 
to Elliot- Yafet dominated regimes when the concentra- 
tion is decreased, as it can be seen in the inset of the 
top panel of Fig. |4] Such crossover is possible when the 
strength of disorder is comparable with the chemical po- 
tential measured with respect to the bottom of the band. 
In this case the spin lifetimes are 3 orders of magnitude 
larger than in the case of hole doping. Note that in the 
electron doped regime a more realistic calculation should 
take into account also the spin splitting of the conduction 
band. 



B. Out-of-plane spin relaxation 

We compute now the out-of-plane spin relaxation rate 
due to an extrinsic or Rashba-likc coupling. In this case 
we have A = —2iXextO'+Sz, which leads to: 



(k) = 4ALt 1 - (a - a') 



2ek 



^ 2 ' 
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FIG. 5: Out-of-plane spin lifetimes as a function of the car- 
rier concentration. In black (dashed) F = 0.001 eV, in blue 
(dotted) F = 0.01 eV, in red F = 0.1 eV. In all the cases 
Aeajt = lO^^Aint- Inset: Spin lifetime for hole concentrations 
where the correction given by Eq. (291. 



The calculation is formally identical to the previous one. 
In the doped regime we have the approximate results: 
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In Fig. [5] the numerical computation of the spin life- 
time as a function of the carrier concentration is shown. 
We take Xext — 10~^Ai„t, which is the correct order 
of magnitude given that this coupling is the result of 
second order processes as we explained in the previous 
section. The spin lifetimes are in this case electron- 
hole symmetric and clearly dominated by the D'yakonov- 
Perel' mechanism. The 1/n behavior is expected from 
the first expression in Eq. (28 1. For /i > A/2 we have 
Tout - hTA'^/i'iTTXl^^aH^nlso for F = 0.02 eV and 
n ~ 10^2 cm~2 we obtain Tout ~ 1 — 2 ns. 

Note that the spin splitting of the valence band is not 
taken into account in this calculation, but its effect is 
relevant since it tends to stabilize the out-of-plane spin 
polarization, in a similar way as an applied magnetic field 
in the out-of-plane direction does. We can take into ac- 
count this effect by cor rectin g the spin relaxation rate for 
hole concentrations asf^^E^ 
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where 2Xint can be interpreted as the Zeeman splitting 
created by an effective magnetic field whose origin is the 
intrinsic SO coupling. Therefore, for F = 0.02 eV and 



10^2 ^ we expect: 



(holes) 
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This correction is taken into account in the calculation 
showed in the inset of Fig. [5] It is interesting to notice 
that our results quantitatively agree with the experimen- 
tal measurements of Mak et al, who have reported spin 
lifetimes exceeding Ins in single layers of MoS2f^ 



IV. DISCUSSION AND CONCLUSIONS 

Our calculations show that the D'yakonov-Perel' mech- 
anism dominates the SO mediated spin relaxation in 
monolayers of M0S2. The method used here is com- 
pletely general, and the results can be extrapolated to 
other dichalcogenides by simply replacing the values of 
the model parameters. Of special interest is the case of 
WS2 monolayers, for which even longer relaxation times 
are expected since it presents an even larger SO coupling, 
of the order of ~ 400 meV.'^SI In general, the in-plane spin 
relaxation is very efficient due to the strong SO interac- 
tion and the lack of inversion symmetry of the system. 
Furthermore it is strongly electron-hole asymmetric due 
to the different spin splitting of valence and conduction 
bands. 

The ah symmetry preserves the out-of-plane spin po- 
larization. However, in real systems this symmetry is 
broken by the presence of a substrate, electric fields, rip- 
ples, etc. Our calculation shows that, for a realistic value 
of the extrinsic or Rashba-like coupling generated by this 
symmetry breaking, the spin lifetimes are of the order of 
nanoseconds. Moreover, the large splitting of the valence 
band due to the spin-obit coupling contributes to sta- 
bilize the out-of-plane spin polarization, as a magnetic 
field in that direction does, and therefore the spin life- 
times in the hole doped regime are expected to be at 
least one order of magnitude larger, in agreement with 
recent experiments.'^"!^ 

It is worth to mention that the symmetry is intrinsi- 
cally broken by the out-of-plane (flexural) phonon modes 



of the M0S2 monolayer, which constitutes an additional 
source of out-o f-plane spin relaxation,!^ as it happens in 
graphene.'SSEl] in the free-standing system the dispersion 
of the flexural acoustic branch is quadratic,!^ as expected 
from symmetry considerations. The bending rigidity of 
the system, which determines the energy of these modes, 
can be estimated assuming a simplified model where the 
M0S2 sheet is described as a plate of certain thickness 6. 
The bending rigidity reads 



Y6^ 



24(l-(j2) 



(31) 



where Y = 0.33 TPa is the Young modulu^^D and cr = 
0.125 is the Poisson ratio. I^Sl If we take as S the inter- 
layer distance, 6 « 6.75 Ajl^Slthen we obtain k « 27 eV. 
Although it is difficult to judge the accuracy of this esti- 
mation due to the dependence, it is reasonable to take 
a bending rigidity bigger than in graphene, even an order 
of magnitude.'^ Therefore, one expects that the higher 
stiffness of M0S2 and the spin splitting of the bands will 
protect the out-of-plane spin polarization. 

In conclusion, we have computed the spin lifetimes of 
monolayer M0S2 within a 2-bands effective model. As- 
suming an extrinsic Rashba-like coupling generated by a 
ah symmetry breaking of the order of Xext ^ 10~^Ai,it, 
we have obtained spin lifetimes of the order of Tout ^ 100 
ns, estimation which is in agreement with recent valley 
population experiments.!^ 
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